Preprint typeset in JHEP style - HYPER VERSION 



RUNHETC-2009-23 



The Imaginary Part of the A/^ = 4 Super- Yang-Mills 
Two-Loop Six-Point MHV Amplitude in 
Multi-Regge Kinematics 



Robert M. Schabinger 

Department of Physics and Astronomy 
Rutgers, The State University of New Jersey 
136 Frelinghuysen Rd. Piscataway, NJ 08854 



Abstract: The precise form of the multi-Regge asymptotics of the two- loop six-point 
MHV amphtude in = 4 Super- Yang-Mills theory has been a subject of recent contro- 
versy. In this paper we utilize the amplitude/ Wilson loop correspondence to obtain precise 
numerical results for the imaginary part of these asymptotics. The region of phase-space 
that we consider is interesting because it allowed Bartels, Lipatov, and Sabio Vera to de- 
termine that the two-loop six-point MHV amplitude is not fixed by the BDS ansatz. They 
proceeded by working in the framework of a high energy effective action, thus side-stepping 
the need for an arduous two-loop calculation. Our numerical results are consistent with 
the predictions of Bartels, Lipatov, and Sabio Vera for the leading-log asymptotics. 



Keywords: ^uper- Yang-Mills Theory, Scattering Amplitudes] . 



Submitted to JHEP 



1. Introduction 

In recent years, tremendous progress has been made towards a more complete under- 
standing of the gluon scattering amphtudes in = 4 Super- Yang-Mills theory (hereafter 
simply N = A). Particularly interesting is the remarkable fact that, in the planar limit, 
it is possible to completely solve the perturbative S-matrix (up to momentum indepen- 
dent pieces) for the scattering of either four gluons or five gluons. Building on the earlier 
work of Bern, Dixon, and Smirnov (BDS) made an all-loop, all-multiplicity proposal 
for the finite part of the maximally-helicity-violating (MHV) amplitudes in = 4. In 
this paper [Q], BDS explicitly demonstrated that their ansatz was valid for the four-point 
amplitude through three loops. Subsequent work demonstrated that the BDS ansatz holds 
for the five-point amplitude through two-loops Q and that the strong coupling form of the 
four-point amplitude has precisely the form predicted by BDS [^. 

In fact, sparked a significant parallel development. Motivated by the fact that the 
strong coupling calculation proceeded by relating the four-point gluon amplitude to a par- 
ticular four-sided light-like Wilson loop, the authors of were able to show that the finite 
part of the four-point light-like Wilson loop at weak coupling matches the finite part of the 
four-gluon scattering amplitude. The focus of was on the four-gluon MHV amplitude, 
but it was shown in that this MHV amplitude/light-like Wilson loop correspondence 
holds for all one-loop MHV amplitudes in = 4. It is clear that an arbitrary n-gluon 
light-like Wilson loop should be conformally invariant in position space. ^ What was not 
at all obvious before the discovery of the amplitude/ Wilson loop correspondence is that 
= 4 scattering amplitudes must be (dual) conformally invariant in momentum space. 

It turns out that this hidden symmetry (referred to hereafter as dual conformal in- 
variance) has non-trivial consequences for the AA = 4 S-matrix. Assuming that the MHV 
amplitude/light-like Wilson loop correspondence holds to all loop orders, the authors of Q 
were able to prove that dual conformal invariance fixes the (non-perturbative) form of all 
the four- and five- point gluon helicity amplitudes (non-MHV amplitudes first enter at the 
six-point level) in AA = 4. Up to trivial factors, they showed that the functional form of 
the (dual) conformal anomaly coincides with that of the BDS ansatz. Subsequently, work 
has been done |8|, ^ that provides some evidence for the assumption made in that the 
MHV amplitude/light-like Wilson loop correspondence holds to all orders in perturbation 
theory. 

The idea is that, due to the fact that non-trivial conformal cross-ratios can first be 
formed at the six-point level, one would naively expect the four- and five-point amplitudes 
to be momentum-independent constants to all orders in perturbation theory. It is well- 
known, however, that gluon loop amplitudes have IR divergences. These IR divergences 
explicitly break the dual conformal symmetry and it is precisely this breaking which allows 
four- and five- gluon loop amplitudes to have non-trivial momentum dependence. In fact, 
the arguments of allowed the authors to predict the precise form that the answer should 

^Strictly speaking, the conformal symmetry is anomalous due to the presence of divergences at the cusps 
in the Wilson loop. If one regulates these divergences and subtracts the conformal anomaly, then the finite 
part of what remains will be conformally invariant. 
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take and they found (up to trivial constants) complete agreement with the BDS ansatz to 
all orders in perturbation theory. 

At this stage, it was unclear whether the appropriate hexagon Wilson loop would still 
be dual to the six-point MHV amplitude at the two-loop level. This question was decisively 
settled in the affirmative by the work of |10| on the scattering amplitude side and |0] on 
the Wilson loop side. Another issue settled by the authors of and [0] was the question 
of whether the BDS ansatz fails at two loops and six points. It had already been pointed 
out by Alday and Maldacena in [12| that the BDS ansatz must fail to describe the analytic 
form of the L-loop n-gluon MHV amplitude for sufficiently large L and n, but it had not 



yet been conclusively proven until the appearance of and |11| that L = 2 and n = 6 
was the simplest possible example of BDS ansatz violation. 

In point of fact, if one was only interested in the question of BDS ansatz violation, it 
would not have been necessary to perform the complete two-loop calculation. It was already 
known to Bartels, Lipatov, and Sabio Vera (BLSV) from computations they performed |13| 
that the BDS ansatz fails at two-loops and six-gluons. It is these computations that we 
will scrutinize in the present paper. For simplicity, we restrict ourselves to the regime of 
multi-Regge kinematics (MRK regime), where all gluons are strongly ordered in rapidity. 
The argument of |13] relied only on the general structure of the multi-Regge limit; it was 
in [0] that the explicit formula for the imaginary part of the analytically continued^ 2^4 
amplitude in the MRK regime was actually presented. 

The goal of the current project is to take the Feynman integral representation of the 
= 4 two-loop six-point MHV amplitude as given^ by ||Tl[| and analytically continue (be- 



fore numerical evaluation) all integrals in the same way that the authors of |14| analytically 
continued their representation of the amplitude. It is then possible to make a direct com- 
parision between exact numerical results and results which should describe the leading-log 
behavior as one goes deep into the MRK regime. The main result of this paper is that, 
when appropriately interpreted, we find solid agreement between our numerical results for 
the imaginary part and the predictions of the approximate formula derived by BLSV. 

In our view this undertaking was necessary because the authors of ||T^, Del Duca, 
Duhr, and Glover (DDG), have expressed doubts as to whether BLSV correctly identified 
the physics of the situation in question. Normally, one expands the amplitude in a Laurent 
series in the dimensional regularization parameter e (defined as d = 4 — 2e) before perform- 
ing any analytical continuations of the amplitude. DDG claim that one can only obtain a 
BLSV-like formula if one goes to the MRK regime and performs all desired analytical con- 
tinuations before expanding the scattering amplitude in e. To the author's knowledge, |16| 
is the only other article that has considered the results of BLSV and DDG.^ In section 8 
of the authors of that work make it clear that they believe the methods of |14| are 

^We postpone a detailed discussion of the phase-space shce under consideration until the next section. 

^As a practical matter, the Feynman integrals appearing in the Wilson loop calculation of [|l^ seem to 
be much easier to evaluate numerically than the Feynman integrals appearing in the scattering amplitude 
calculation of p^. 

^In the first version of this article remarks were made about that have been removed in this version. 
The reason for this change is that, in going from v3 of to the final published version, the discussion of 
BLSV and DDG changed significantly. We thank Howard Schnitzer for bringing this to our attention. 
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correct. Their paper, however, focuses on somewhat different issues and does not address 
the question quantitatively. Our goal is to demonstrate conclusively that there are no 
problems with the work of BLSV (at least for the 2 — > 4 amplitude that we have consider 
in detail) and it is hoped that our results will put to rest the controversy surrounding the 
results presented in [p!^]. 

This article is organized as follows. In section ^ we establish our notation and describe 
in more detail the problem which we have solved. In section ^ we present the main results 
of this paper and compare them to those of BLSV. In section Q we present our conclusions. 



explain how our work sheds light on some claims made in reference |15|, and discuss some 
open problems appropriate for future work. In Appendix]^ we review six-body kinematics 
and the relation between scattering amplitude and Wilson loop notation. Finally, in Ap- 
pendix we give some formulas useful for analytical continuation and an explicit example 
of one of the more intricate calculations we performed. 



2. Preliminaries 



As mentioned in the introduction, we do not deal with the two-loop six-gluon MHV ampli- 
tude directly. Instead, we use the two-loop light-like hexagon Wilson loop which is known 
to be equal to the amplitude of interest |ll|] in a sense that will be made clear below. 
Recall that, for a given multiplicity, there is effectively only one MHV amplitude in = 4, 
due to the power of supersymmetric Ward identities (see reference [0] for a review of these 
ideas) . This fact suggests that it is most natural to consider MHV loop amplitudes divided 
by the appropriate Parke- Taylor |jl8| tree amplitude. The way in which one formally ex- 
presses the correspondence between the finite parts of 2-loop six-gluon helicity amplitudes 
and the 2-loop light-like Wilson loop is 



In 



' <=^(l-,2-,3+,4+,5+,6+) ' 
At'--(l-,2-,3+,4+,5+,6+) 



finite; 0{a'2) 



In {A^>^ 



ht-likc WL^ 



(2.1) 



finite; 0{a2) 



where we have chosen the (1 ,2 , 3^, 4+, 5^, 6"*") helicity configuration for sake of definite- 
=>5 What we really mean by the above equivalence is that, if the logs on either side are 



ness 



expanded to 0{a^) (the 't Hooft coupling, o, is defined below), the finite parts of the 0(0^) 
coefficients should match up to an additive constant independent of the dynamics. The 
additive constant which turns the above equivalence into an equality was given in jl9[|.^ 

Before we try to make sense of Eq. (|2.l| ), let us first remind ourselves of the overall 
normalization used for the numerical evaluation of the finite part of the hexagon Wilson 



loop. In 1 11], the exact color structure was computed and only at the end of the calculation 
did the authors take the Nc ^ oo limit. They were then able to pull an overall factor of 
CpNcig'^ /4:Tt'^)'^ out of the finite part of all Feynman diagrams. Finally, they made the 



^Note that this does not imply 



(l-,2-,3^ 



an equivalence between 



finite parts of 
See references 118(1 



the 0{a^ 

,4+,5+,6+) divided by the tree amplitude and (^^ight-iike wl^ 
and for details. 

®It should also be pointed out that the 0(e^^) pole terms do not exactly match. The physics of this 
mismatch was properly understood in |2l[]. 



-3- 



replacement CFNc{g'^ /^ir'^)'^ — > 2a^. This replacement is that which is required to obtain 
numerical results consistent with the definition used by BLSV for the 't Hooft coupling: 



87r2 



(2.2) 



As mentioned above, proceeding naively results in a mismatch between the two-loop 
hexagon Wilson loop and the two-loop six-point MHV amplitude. To understand this 
problem, we need a little more information about the form of the expected answer. On the 
amplitude side, the answer is constrained by dual conformal invariance to look like 



In 



'<=4(1- 



,2-,3+,4+,5+,6+) 



^t'--(l-,2-,3+,4+,5+,6+) 



77BDS 



(2.3) 



finite; Oia^) 



where -^^2*^ is the finite part predicted by BDS in Q (we use the conventions of l]!^] for this 
function) and the Ui are conformal cross-ratios built out of the standard cyclicly symmetric 
Mandelstam invariants (see Appendix |A| for a discussion of the essential points): 



ui 



S1S4 



U2 



S2S5 

tlt2 



U3 



t2h 



(2.4) 



R&;2{ui-,U2-,Uz) is simply called the (two-loop six-point) remainder function. By con- 
struction, Rq-^2{uitU2tU3) ^ in any collinear limit. Operationally, this means that 
Rq;2{ui,U2,uz) ^ as Uj ^ 1; Uj+i, Ui+2 ^ (i mod 3). 

However, if one uses the dictionary between Mandelstam invariants in the scattering 
amplitude picture and non-zero lengths in the Wilson loop picture (Appendix |A|) , one finds 
that, on the Wilson loop side, the answer looks like H 



In {A, 



light- like WL|^ 



^6^2°^+ (^6;2(nl,n2,n3) + — 



finite; 0(0^) 



vr 



(2.5) 



Therefore, before we can utilize the Wilson loop picture calculations of 1 11], we have to ad- 
just their formulas by subtracting an overall factor of vr'^/S on the right-hand side. This gu- 
rantees that the Wilson loop version of R&-2{u\,U2-,U3) has the correct (vanishing) collinear 
limit. 

We now describe precisely what the MRK regime is in the context of this paper. 
Typically, when one performs a loop-level gauge theory calculation one takes all invariants 
to be negative and only analytically continues the result back into the physical region after 
all Feynman integrals have been evaluated''. The region of phase-space with all invariants 
negative (hereafter referred to as the Euclidean region) is special because all Feynman 
integrals evaluate to real numbers. We will eventually be interested in a specific analytical 
continuation, but for now we discuss the Euclidean MRK regime. 



^Four-body scattering is a good example. Physically, one is typically interested in a 2 — * 2 process with 
s > and t < 0. On the other hand, for the evaluation of Feynman integrals, it is convenient to take all 
particles outgoing and both s and t negative. 
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The Euclidean MRK regime is one in which all of the Mandelstam invariants are 
space-like and their absolute values satisfy the following ordering relations 



Sl » S5,Si,S3 » S6,t2,S2 , (2.6) 

where we have suppressed the absolute value bars for the sake of simplicity. Any phase- 
space slice that is conformally equivalent to one satisfying (|2.6D is suitable for studying 
the remainder function in the MRK regime. In fact, for this purpose one can define the 
MRK regime completely in terms of ui, U2, and u^. In this regime all of the particles are 
collinear (either left or right movers) and, consequently, two of the conformal cross-ratios 
approach zero (from above) and one of them approaches unity (from below). Which cross- 
ratio approaches unity depends on how one labels the external momenta. The particular 
labeling scheme that we have made in this work coincides with that of and gives ui ^ 1 
and U2, M3 ^ 0. In any case, the fact that the remainder function is unchanged by cyclic 



permutations of the cross-ratios |11] ensures that all choices lead to the same physics. We 



refer the interested reader to the introduction sections of [Q, and [^] for a much 
more detailed discussion. As expected, we find that the remainder function in the Euclidean 
MRK regime is approaching zero. Due to the fact that the remainder function approaches 
zero in the Euclidean MRK regime, it is not suitable for detecting BDS ansatz violation at 
the two-loop six-point level. 

Instead we consider a particular analytical continuation of the Euclidean MRK regime 
that preserves the conformal cross-ratios. The analytic continuation is one where both 
Mandelstam invariants in the numerator of ui, the cross-ratio approaching unity, are al- 
lowed to change sign. The nice thing about such an analytical continuation is that one 
can set up the problem in the Euclidean MRK regime and then analytically continue by 
simply making the replacements —si si and — S4 — > S4. BLSV claim that, even though 
R6]2{ui,U2,U3) ^ in the Euclidean MRK regime, the remainder function develops a 
non-vanishing imaginary part (that grows logarithmically) if one performs the analytical 
continuation described above, thus providing a slick way to check that the BDS ansatz fails 
at the two-loop six-point level. We evaluate this claim quantitatively in the next section. 



3. The Imaginary Part 

In principle there is no need to restrict oneself to the calculation of the imaginary part, but 



we do so both because this is what BLSV presented in |14] and because the imaginary part 
turns out to be much easier to control numerically. To make the calculation as simple as 
possible, we attempted to parametrize the non- vanishing distances in a way that makes the 
number of Feynman integrals that need to be analytically continued as small as possible. 
Unfortuately, some complication is unavoidable, since several of the higher-dimensional 
Feynman integrals have to be processed further in the Euclidean region (beyond what was 
done in the Appendix of |jll|]) to put them in a form suitable for analytical continuation. 
Specifically, we found that, at least for numerical evaluation in Mathematica, all Feynman 
diagrams had to be expressed in terms of integrals which were two-fold or less. Before 
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continuing, we first made sure that all of the contributing integrals were correctly entered 
into Mathematica by reproducing Table 1 of |11]. 

The primary phase-space slice that we worked on can be parametrized in terms of the 
three positive real numbers a, b, and c and is given by 

2_ 2_ 2_i,2_ 2_ 2_i, 

3^13 — c, — a, X35 — 0, — c, X]^5 — a, — 

= —1; X25 = —1; x\q = —1 . (3.1) 

On this slice we have the cross-ratios ui = (? , U2 = 0? , and 1*3 = 6^, where 



c = a + 6 - 1 + 2VI -a-b + ab (3.2) 

is fixed by the Gram determinant constraint (see Appendix^ for details). 

Analytical continuation to the region of interest is performed by replacing — c with c 
and then carefully following the rules outlined in Appendix^. We found that, apart from 
several trivial factors similar to those that appear in the BDS ansatz, we had to continue 
84 integrals away from the Euclidean region to extract the complete imaginary part of 
-R6;2(c^) fl^) &^) on this particular branch. 

Carrying out the calculation is quite arduous, due to the fact that Mathematica is 
not equipped to implement the rules described in Appendix ^ automatically. Fortunately, 
conformal invariance can be used as a stringent cross-check on the final results. One has 
to look at a different kinematic point that has the same conformal cross-ratios and make 
sure that Rq-2{c'^ , , b'^) is defined on the same branch after the analytical continuation 
(— c — > c) has been performed. We generated a secondary test configuration along these 
lines given by 

2_ 2_ 2_26 2_o 2_ 2 _ nj 

^13 ~ "'^i ^24 — ^35 ~ ""3"' "^^6 ~ -^15 ~ -^26 ~ ^" 

^14 ~ ~2' ~ ~3' ^36 (3-3) 

It is easily checked that this configuration preserves the cross-ratios. Analytically continu- 
ing this configuration in the manner described above, we confirmed that the result for the 



imaginary part of i?6;2(c^, o^, ^^) on (3.3) is equal to the imaginary part of Rq-2{c'^ , a'^ , b"^^ 



on ( p. ID to within the available numerical precision. 

Once we established confidence that our Mathematica script for the primary configu- 
ration was correct, we extrapolated our results deeper into the MRK regime and compared 



them to the explicit formula given by BLSV |14| 



lm\Re,2iui,U2,U3)} = J In (1 - ^xi) In (^^) In (^^) + • • • (3.4) 



in terms of dual conformal cross-ratios (we've pulled out two factors of a to make contact 
with the conventions of |jll[). The dots represent sub- leading- log terms not determined 



by the analysis in |14| 



At this stage the skeptical reader may be wondering how we intend to make a mean- 
ingful statement about the BLSV function using numerical analysis when the sub-leading 
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corrections to it are completely unknown. There are really two issues that need to be ad- 
dressed: Does the BLSV function capture the leading-log behavior of the imaginary part? 
If so, how much do the sub-leading corrections matter? In general, it is not guranteed that 
numerical analysis can be used to address these questions. Fortunately, the problem at 
hand has a very special structure which allows us to make progress. 

Before discussing our results, we have to explain how we parametrized the problem. 
We were motivated by the structure of the Gram determinant constraint to think about 
expanding the expression for mi(m2,^*3) in {u2,U3) about (0,0). This allows us to Taylor 
expand Eq. (|3.2|) 



1 — Ui 



(3.5) 



and drop the higher-order terms. In this approximation, we can rewrite the BLSV formula 
in terms of just two variables, which we choose to be U2 and ^ = U2/U3: 



Im|i?6;2(u2,0} = "^^"^ ( ^2 



1 



1 



In 



*)■ 



In 



+ 



(3.6) 



Note that Eq. ( ^.61) depends on U2 through the first logarithm only. This means that on 
lines of constant ^ the BLSV function looks like 



Aln{u2) + B + 



(3.7) 



with A and B given by 




(3.8) 
(3.9) 



plus the unknown subleading-log contributions which one might think could depend on U2 
and U3 in a complicated way, even for U2 and U3 very close to zero. 

In fact, it turns out that the sub-leading terms depend only on ^ to leading order. If the 
sub-leading logs depend only on (a fact that we justify a posteriori), then grouping -B(^) 
with the A(S^) ln(n2) term in the BLSV function is an arbitrary choice; for the purposes of 
this analysis we could equally well group it with the unknown corrections because it too 
depends only on ^. Moreover, the assumption that the unknown contributions depend only 
^ suggests a simple experiment to see whether the BLSV formula actually does capture the 
leading-log behavior of the imaginary part. For a given ^, we take nine data points evenly 
spaced on a log scale (6.9444 x 10"'^ < ii2 ^ 0.0016), differentiate^ with respect to ln(n2) 



^To differentiate our data with respect to ln(u2) we took two data points very close together and 



computed U2 



Au2 



for the pair. We then experimented with different seperations in U2 to 



make sure that the results were stable. 
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and then plot the results (Figure ||) against Eq. ( |3.8D as a function of ln(u2). The range of 
ln(«2) plotted was fixed for us by the breakdown of our numerics for very small U2 and the 
breakdown of our approximation scheme (Eq. (^]^)) for relatively large U2- We performed 
the analysis for four evenly spaced values of the ratio ^ on the unit interval. Recall that 
the imaginary part is symmetric under the interchange of U2 and u^. This implies that we 
are free to fix U3 > U2 without loss of generality. In this case, (0, 1) is the natural domain 
for ^.^ 

We find that the numbers extracted from our data match those from Eq. ( |3.8D very 
well indeed. The fact that this works at all justifies our assumption that the unknown 
sub-leading-log contributions depend only on ^ to first order in the small quantities U2 and 
U3. Furthermore, it shows that our exact results are consistent with a non-trivial prediction 
of BLSV as a function of ^. It should also be stressed that, although, we produced values 
of ui by solving the Gram determinant constraint in the Euclidean MRK regime for our 
primary kinematic configuration, this was not a necessary step and we could have obtained 
qualitatively identical results without the Gram determinant constraint. After all, our 
secondary test configuration did not satisfy the Gram determinant constraint and yet, by 
virtue of the conformal invariance, it matched a configuration that did satisfy the constraint. 
For example, if we scaled 1 — ui to zero at exactly the same rate that U2 and were scaled 
to zero we would have been able to produce plots qualitatively equivalent to those displayed 
in Figure |l|. 
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Figure 1: The derivative of our data with respect to ln(M2) (black dots) and the coefficient A of 
Eq. ( ^3.81) (red dashed) are plotted as a function of ln(M2). From bottom to top, results for U2/M3 (C) 
ratios 0.2, 0.35, 0.6, and 0.8 are shown. Note that the lowest-order approximation used for 1 — ui 
(Eq. ( |3.5| )) starts to break down for ln(u2) greater than about -8. 

Of course, there is a piece of Eq. (|3.6D (the function B{^) given by Eq. ( ^.9D ) that 

^At u-j, = U2 there is a branch point and for sufficiently large U3 (sufficiently small ^) we will necessarily 
leave the MRK regime if we consider contours of constant ^ over the range of ln(M2) plotted in Figure H. 
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cannot be detected by taking derivatives with respect to ln(it2)- We expect that it is 



present by virtue of the fact that the analytical continuation performed in |14] was in the 
invariants si and S4 sitting inside of ui, but we have no way of nailing it down. Although 
we can't say much about -B(^) or the unknown ^ dependent corrections, it would be nice 
to get a qualitative feel for how close the BLSV function comes to describing the complete 
imaginary part. Of course, the phase-space slice considered above does not allow one to 
do this because two of the three logarithms in the formula were nearly invariant on it. We 
can do better by varying ^ keeping U2 fixed. We plotted our data for U2 = 1.5625 x 10~® 
(an arbitrarily chosen point that was found to be numerically stable and solidly within the 
regime of validity of the small (u2, U3) approximation) and varied M3 from 1.0609 U2 to 36 
U2. These results are compared to the BLSV function and A(^)ln(u2) from Eq. (|3.7| ) in 
Figure |2|. We see that the general shape of the curves is predicted pretty well by BLSV 
but that the sub-leading logs become progressively more important as the branch point at 
^ = 1 is approached. The fact that A{S^)ln{u2) (the blue dot-dashed curve) fits the data 
plotted in Figure || better than the actual BLSV formula (red dashed curve) suggests that 
the true leading-log behavior is described by A(S^)ln{u2) on the phase-space slice under 
consideration (Eq. ( p.lD ) and that it is more natural to group -B(^) with the undetermined 
sub-leading-log contributions. 
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Figure 2: Our data (black dots), the BLSV function (red dashed), and A{^) ln(u2) from Eq. (3.7) 
(blue dot-dashed) are plotted as a function of U2/W3 (^) for U2 — 1.5625 x 10^^. Presumably, the sig- 
nificant deviations present for large ^ (for both curves) are due to the presence of the undetermined 
sub- leading- log contributions to the imaginary part. 



4. Discussion and Conclusions 

In this work we studied the claim made by BLSV that a particular analytical continuation 
of the Euclidean MRK regime allows one to see that the BDS ansatz fails for the two-loop 
six-gluon MHV amplitude in = 4. We find that our precise numerical results for the 
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imaginary part of the amplitude are completely consistent with the formula derived by 
BLSV in |14|. Our work provides strong evidence that the interpretation of the physics 
given in the work of BLSV is correct, at least for the 2^4 process that we considered in 
detail. Our work also bears on a couple of issues raised by DDG. 

In |jl5|, DDG claim that it is not possible to see the imaginary part reported by BLSV 
unless one analytically continues away from the Euclidean MRK regime before expanding 
the amplitude in e. We encountered no such difficulty. In other words, we have no reason 
to believe that there is any problem with the conventional approach to perturbative gauge 
theory calculations in this context. Furthermore, DDG remark that the kinematical regime 
chosen by BLSV is unphysical. It is not clear to us what DDG mean by this statement. 
In any case, we show in Appendix ^ that, if desired, it is possible to solve the Gram 
determinant constraint after analytical continuation (— c c) on the phase-space slice 
defined by Eq. Q. 



Finally, we highlight a point of potential confusion, also discussed in [16|. One might 
worry that the analytical continuation described in Section ^ is ill-defined because the so- 
lutions to the Gram determinant constraint on a phase-space slice like the one defined in 



Eq. (3J) are not invariant under — c — > c (compare A. 4 and |A.5| to A. 6 ). Our analysis has 
explicitly shown that a well-defined imaginary part exists regardless of whether the Gram 
determinant constraint is satisfied. This is most clearly seen by recalling that we were able 
to obtain results using ( |3.3D identical to those found by using ( |3.1| ), even though the kine- 
matic configuration defined by (|3.3| ) did not satisfy the Gram determinant constraint before 
or after the analytic continuations were performed. Several authors |]l^, |l^] have already 
concluded that the (dual) conformal invariance of the remainder function is independent 
of the Gram determinant constraint. Although their calculations were performed in the 
Euclidean region, our work has shown explicitly that the imaginary part of the remainder 
function behaves in a completely analogous way. 

Our main goal in this endeavor was to bring some clarity to the existing literature. One 
deficiency of the present paper is that the analysis we have presented is purely numerical 
and is therefore somewhat limited in scope. It may prove useful to find an analytical 
formula for the imaginary part we have studied, because this result could yield useful 
information about the functional form of the remainder function for general kinematics. It 
is also quite likely that this line of attack will be simpler than tackling the general problem 
directly because the space of transcendentality three functions is smaller than the space of 
transcendentality four functions. 
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A. Six-Body Kinematics and the Amplitude/ Wilson Loop Correspondence 

In this appendix, we review some facts about six-body kinematics and estabhsh a dictionary 
that makes the mapping between amphtudes and Wilson loops concrete. 

Although there are only eight independent Mandelstam invariants in a six-body scat- 
tering process [^3| , it is more convenient for our purposes to use a set of nine variables that 
arise naturally in the propagator denominators of Feynman diagrams contributing to the 
scattering process. Taking all external momenta to be outgoing, these variables are defined 
by the relations 

Si = {h + ki+if and tj = {kj + kj+i + kj+^f , (A.l) 

where the index i is taken mod six and then index j is taken mod three. This basis of 
Mandelstam invariants has the nice feature that it forms a cyclicly symmetric set. 

There are several references in the body of this paper to the "Gram determinant 
constraint" . This constraint is what makes only eight of the nine variables defined above 
independent quantities. The constraint is simply expressed pS] as 



Dei' {2ki ■ kj) = , (A.2) 

where the prime signifies that one of the six four-vectors in the problem is to be omitted 
in taking the determinant. To be concrete we omit k^. In this case, we have 

{2ki ■ kj)k, = (A.3) 

/ Si ti - Si - S2 S2 + - ti - t2 t2- S5- S6 \ 

Si S2 t2- S2- S3 S3 + Se - t2 - h 

h- Si - S2 S2 S3 t3- S3 - Si 

S2 + S5 - ti - t2 t2- S2- S3 S3 S4 

V ^2 - S5 - S6 S3 + Sq - t2 - t3 t3- S3 - S4 S4 / 

There are actually two ways that one could set Det(2A;j-/cj)fcg equal to zero for the kinematics 



defined by Eq. (3.1) and still be in the Euclidean MRK regime. Although we chose the 
solution 

c = a + b-l + 2VI - a-b + ab (A.4) 
for the analysis performed in this paper, we equally well could have chosen 



Vl-a^-b^ + - ab (A.5) 



because both ( |A.4| ) and ( A.5 ) force c to approach 1 from below when a and b are taken 
to zero from above. It is also worth noting that, if desired, one could solve the Gram 
determinant constraint in the MRK regime after analytical continuation (— c — > c). In this 
case, there is a unique consistent choice given by 



\/l - a2 - 62 + ^252 + , (A.6) 
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We intentionally worked with kinematics where the Gram determinant constraint was 
satisfied in the Euclidean MRK regime before analytical continuation (c(a, b) defined by 
Eq. ( [A.4D ) because we wanted to see explicitly that it does not affect the analysis if the 
Gram determinant constraint fails to be satisfied after the replacement — c — > c has been 
made. 

Finally, the mapping between distances in the Wilson loop picture and the canonical 
basis of Mandelstam variables defined above is expressed as 



■'^13 


<— > 


Si 






S2 


•'^35 




■S3 


^2 


•5— > 


S4 


2^15 




•S5 


^2 


<— > 




^2 


<— > 


h 


^2 


<— > 




™2 
■^36 


<— > 





B. Analytical Continuation 

In this appendix, we discuss the essentials of analytical continuation as it pertains to the 
results presented in this paper. As an example, we give a complete derivation of the 
imaginary part for the integrals associated to topology (j) in [0]. All of the formulae 
below are easily derived from the results given in Appendix C of |25]. Here we simply give 
a representative sample of the results needed for the particular problem studied in this 
paper. In what follows, s and t should be thought of as negative real numbers such that 
St < 1. 

Inf— ^ =lnf— ^ -iTT (B.l) 



tj \-t^ 

ln(j) =ln(5£) (B.3) 

LMl-^) .^_l,.,(,,Z£)_L,(^)+i.,„(l + Zf) ,B.4) 

U2{l-j) ^uJl-^\ (B.6) 

LijIl-sA =Li,(l-(-s)(-t)) + 2iir 111(1- (-»)(-*)) (B.7) 
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= Li2( 








= Li2f 



^ 2i7rln(l-— ) (Bi 



^ + 2tt^ - Iniln ( ^- — r-l) (B. 



{-s){-t)J \{-s){-t) 



The starting point for topology (j) is Eq. (A. 23) of |11]. Though we technically have 
a Wilson loop integral, we express everything in terms of Sj and tj using the dictionary of 
Appendix |^ because it makes the formulas easier to read. It is also convenient to define 

N{si,S4, ti,t3,Vi,V2) = ts - (^3 - 54)^^1 - {ts " Si)v2 - (si + S4 - ti - t3)viV2. (B.IO) 

After a little reorganization, we pull out two factors of the 't Hooft coupling and sum over 
distinguishable cyclic permutations to obtain 



, 1 , f N{si,S4^,ti,t3,Vi,V2) 

dv2 7—7^ ■ In 



1 1 f N{si,S4,ti,ts,Vi,V2)\ , 

X 7 In ; ^ + 

-Jf-l^+V, V h{l-V2)+S,V2 J 

1 In ( ^('^2,g5,^2,tl,^l,f2) \ ^ 



1 / 2 +^ Qo 



,3+|_y3-t. +^2 V t2(l-^^2)+S3^2 . 

On the slice defined by Eq. (p. ID, this becomes 

/ Af(-c,-c,-l,-l,i>i,i>2) \ ( N(-c-c-l-l,v-i,V2) 

1 /"l I C(l — V\)+Vl j \ 1-V2+CV2 

(l-2.,)(l-2..) ^ + 

, / N{-a,-a-l-l,V'i_,V2) \ , / Af(-a,-a,-l,-l,-i;i,ii2) 
1 /■! I a(l-,;i)+i;i I l-U2+at'2 

Jo (1 - 2z;i)(l - 2z;2) 

, / jV(-fe,-b-l -l,Di,i;2) A i„ ( N{-b,^b -1,-1 m,V2) 
1 /-l I 6(l-i.i)+t;i I \ l-V2+bv2 



Jo (1 - 2t;i)(l - 2t;2) 

Clearly, only the first of these terms develops an imaginary part under — c — > c. After 
dropping the other terms and making this substitution, we find that we are faced with 
analytically continuing 



. In f m££rzhz^iml\ In ( N{c,c-i,-i,v,,v2) \ 

fl /"i 1 -C{l-Vl)+Vl / 1 1-V2—CV2 / 

Jo Jo (l-27;i)(l-27;2) 



(B.13) 



-13- 



on the unit square in the vi — V2 plane. A glance at Eqs. (|B.1|) , ( [B.2D , and ( |B.3|) tells us 
that we must isolate the regions of this parameter space on which either the numerator or 
denominator (but not both simultaneously) of either logarithm change sign. The conditions 
for the denominators to change sign are simple because they are expressed in terms of 
constants. For the denominator of the first logarithm the condition is 

n < (B.14) 

and for the denominator of the second logarithm the condition is 

The condition for the numerator is more complicated. We choose to express the conditions 
as a function of V2 and they are then given by the inequalities 

1 - V2{l + c) r 

Vl > — — II V2 < 



(l + c)(l-2^;2) 1 + c 

l-?J2(l + c) -r. ^ 1 ^-D^r\ 

' - (l + c)(l-2w2) '-1 + C ^ ^ 

It isn't possible for the numerators to change sign if V2 is between c/(l + c) and 1/(1 + c). 

We now systematically investigate how the numerators and denominators of the log- 
arithms in ( |B.13| ) change sign as one moves between different domains. Let us begin 
by parametrizing the problem. In V2 there are three regimes: < V2 < c/(l + c), 
c/(l + c) < f2 < 1/(1 + c), and 1/(1 + c) < t;2 < 1- In the first, we find three differ- 
ent regions: < vi < c/(l + c), c/(l + c) < wi < (1 - ^2(1 + c))/(l + c)/(l - 2v2), and 
(1 — f2(l + c))/(l + c)/(l — 2^2) < Vl <1. In the second, we find just two: < vi < c/(l + c) 
and c/(l + c) < Vl < 1. Finally, in the third regime, we again find three separate regions: 
< < (1 - ^2(1 + c))/(l + c)/(l -2v2),a- V2il + c))/(l + c)/(l - 2v2) < vi < c/(l + c), 
and c/(l + c) < i;i < 1. 

We now describe what happens in each of these eight regions. If < f 2 < c/(l + c) and 



< < c/(l + c), then the denominator of the first logarithm in ( B.13 ) changes sign. If 
< V2 < c/(l+c) and c/(l+c) < vi < (l—?;2(l+c))/(l+c)/(l— 2^2), then nothing happens. 
If < f2 < c/(l + c) and (1 — ^2(1 + c))/(l + c)/(l — 2^2) < f 1 < 1, then both numerators 
change sign. If c/ (1 + c) < ^2 < 1/(1 + c) and < f 1 < c/(l + c), then the denominator of 
the first logarithm changes sign. If c/(l + c) < V2 < 1/(1 + c) and c/ {1 + c) < vi < 1, then 
nothing happens. If l/(l + c) < i;2 < 1 and < wi < (1 - ^2(1 + c))/(l + c)/(l - 2?;2 ), then 
all numerators and denominators change sign, effectively the same as if nothing at all had 
happened. If 1/(1 + c) < ^2 < 1 and (1 - V2{1 + c))/(l + c)/(l - 2^2) < vi < c/(l + c), 
then both denominators change sign. Finally, if 1/(1 + c) < V2 < 1 and c/ (1 + c) < f 1 < 1, 
then the denominator of the second logarithm changes sign. 



Putting all of this information together, one can see that ( B.13 ) is easily expressed in 
terms of four double integrals: 

1 c I In ( ■'^('^'C ] 4- i-jr] In ( ^(g.c.-i,-i,-»i,-»2) 

T+S f-rh \ \ c(l-vi)-vi j ^ j \ 1-V2-CV2 

dv2 I dv-[ 
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(l-27;i)(l-2?;2) 
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1 + c 



+ / dvo dv 



1 \ '''\ -c(l-Vl) + Vl ) M \ 1-V2-CV2 



+ / dv2 , , dv 



l + c (l + c){l-2i;2) 

.1 .1 i°( "-:i7i-;:srO ('"( "''-.^i;:^r' )+" 

(i-2;.)(i-2..) — <'''^* 

l + c l + c 

The imaginary part of the above expression is easily extracted and is equal to 

^ c In ' ^(^'g'-l'-^'^i'^a) 

in / dv2 / dv 
Jo Jo 



1~V2—CV2 

1 (l-2t;i)(l-2t;2) 



I -iV(c.c -l,-Mi.^2) \ ^ [ -Nic c,-l,-l,v„V2) 

i I — C(l — Dlj+fl / I 1 — V2—CV2 



_L^^2li+£) (1 - 27;i)(l - 2^72) 

(l+c){l-2u2) 

Af(c,c-1 -l,i;i,i>2) \ _^ 1^ Af(c,c,-l,-l,t)i,t;2) 



-1 rjf^ \ c(l-vi)-vi I ^ \ -1+V2+CV2 

-in I dv2 / dvi 



i-^2(i+c) (1 - 2t;i)(l - 2u2) 



l+c (l+c){l-2u2) 

1„ / A^(e,c,-l,-l,fi,D2) 

+in / dua / dvi—^ . (B.18) 

' ■ ' (l-2z;i)(l-2t;2) ^ ' 



1 

l + c " l + c 



This particular topology works well as an example because it is compact enough to be 
presented in detail and, nevertheless, it is one of the most intricate continuations that one 
must perform to reproduce the results given in this paper (working on the slice defined by 
Eq. Q). 
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